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We report on the nonequilibrium quasiparticle dynamics in BaFe2As2 on both the hole doped 
(Bai_ x K x Fe2As2) and electron doped (BaFe2- y Co y As2) sides of the phase diagram using ultrafast 
pump-probe spectroscopy. Below T c , measurements conducted at low photoinjected quasiparticle 
densities in the optimally and overdoped Bai_ x K x Fe2As2 samples reveal two distinct relaxation 
processes: a fast component whose decay rate increases linearly with excitation density and a slow 
component with an excitation density independent decay rate. We argue that these two processes 
reflect the recombination of quasiparticles in the two hole bands through intraband and interband 
processes. We also find that the thermal recombination rate of quasiparticles increases quadratically 
with temperature in these samples. The temperature and excitation density dependence of the 
decays indicates fully gapped hole bands and nodal or very anisotropic electron bands. At higher 
excitation densities and lower hole dopings, the dependence of the dynamics on quasiparticle density 
disappears as the data are more readily understood in terms of a model which accounts for the 
quasiequilibrium temperature attained by the sample. In the BaFe2- y Co y As2 samples, dependence 
of the recombination rate on quasiparticle density at low dopings (i.e., y = 0.12) is suppressed 
upon submergence of the inner hole band and quasiparticle relaxation occurs in a slow, density 
independent manner. 

PACS numbers: 74.25. Gz, 78.47.-fp 



I. INTRODUCTION 

The discovery of the iron-pnictide superconductor^^! 
has opened a new chapter in the study of high- 
temperature superconductivity. Although they display a 
lower degree of electronic correlation than the cuprates^, 
there are a number of similarities between the two mate- 
rial classes. Perhaps the most striking of these is the close 
correspondence between their phase diagrams. In both 
classes, superconductivity is achieved by the doping of 
either electrons or holes into a "parent" material which 
displays magnetic ordering at lower temperatures, i.e., 
spin density wave (SDW) ordering in the pnictides^! or 
antifcrromagnetic ordering in the cuprates. As this static 
magnetic ordering is suppressed with doping, supercon- 
ductivity emerges, increasing in transition temperature 
T c to a maximum at optimal doping before receding with 
the addition of further carriers. 

With these shared characteristics come shared myster- 
ies, including the natures of the pseudogap state, of the 
coexistence between magnetic and superconducting order 
(i.e., homogeneous vs. heterogeneous), and, arguably the 
most fundamental, the binding boson itself. The pres- 
ence of a resonance in neutron scattering spectra below 
T c in the pnictides^! and the suppression of magnetic 
ordering upon the addition of carriers suggests that an- 
tiferromagnetic spin fluctuations mediate pairing, imply- 
ing that the superconducting order parameter A switches 
sign in the Brillouin zone. Therefore, definitive determi- 
nation of the gap structure and symmetry would serve 
as a significant step towards understanding the mecha- 



nism of superconductivity. In the cuprates, experiments 
clearly show d-wave pairing^. The picture in the pnic- 
tides is a lot murkier, in part due to their multiband char- 
acter and the fact that interband interactions likely play 
a role in superconductivity^^. This complexity may be 
at the origin of the seemingly conflicting sets of experi- 
mental evidence concerning the symmetry of A. Angle 
Resolved Photoemission (ARPES) studies have consis- 
tently revealed a fully gapped Fermi surface^E^ w hile 
NMlP^ and magnetic penetration deptlP^ measurements 
suggest the presence of nodes. Confronted with these 
data, theorists have proposed a diverse collection of or- 
der parameters, including an s± symmetry with nodes^, 
d xy symmetry^, and the prevailing nodeless s± order 
parameter with interband impurity scattering^. 

In this study, we use optical pump-probe spectroscopy 
to investigate the gap symmetry and interband interac- 
tions in the Ba-122 system of the iron pnictides. Pump- 
probe spectroscopy is a powerful technique in which ab- 
sorption of an ultrashort "pump" pulse results in the in- 
jection of a transient population density n of nonequi- 
librium quasiparticles^. A second "probe" pulse then 
records the return of the system to equilibrium through 
measurement of the time-resolved change in reflectivity 
(AR(t)/R), assumed linearly proportional to n. These 
measurements have been successfully employed in the 
cuprates to examine, e.g., gap s ymme tr y! 20 ! 21 ! and the 
strength of electron-boson coupling} 22 ! 23 !. i n the pnictides, 
pump-probe spectroscopy has probed the existence of a 
pseudogap state^, the competition between SDW and 
SC ordering^, coherent lattice vibrations to rule out 
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their role in superconductivity^, and the role of inter- 
band interactions and gap symmetry in nonequilibrium 
quasiparticle relaxation^. As prior ultrafast experiments 
have focused on one side of the phase diagram at a time, 
there is a need for a systematic study of the pnictidcs as 
a function of doping across the phase diagram. 

Building upon our prior worlPH, we study both 
the hole-doped Bai_ :E K :r Fe2As2 and the electron-doped 
BaFe 2 _j,Co a As2 series of the 122 pnictides. In the opti- 
mal to overdoped hole doped samples (x = 0.4,0.5,0.6), 
the signal (AR(t)/R) is composed of two distinct fea- 
tures. The first is a fast component whose decay 
rate de pends o n excitation fluence <E> which, aided by 
ARPE gi2ii4i28| measurements and LDA calculations^), 
we argue arises from the inner hole bands. The second 
is a slow, fluence independent decay we attribute to the 
outer hole band. An ana lysis of the fast component based 
on the Rothwarf- Taylor- 20 * 30 ! coupled differential equa- 
tions reveals a T 2 dependence of the thermal population 
of quasiparticles on temperature, suggesting the presence 
of nodes on the Fermi Surface. At the same time, the 
number of photoinduced quasiparticles that contribute 
to the overall signal increases linearly with laser fluence, 
reflecting the fully gapped character of these bands. We 
thus argue that these observations are consistent with 
fully gapped hole pockets at zone center and nodal or 
highly anisotropic electron pockets at the zone bound- 
ary. 

In the normal state, we observe oscillations in the re- 
flectivity transients due to coherent acoustic phonons. 
These thermally driven acoustic waves are suppressed be- 
low T c due to a significant decrease in the thermal expan- 
sion coefficient below the transitiorP^H _ except at the 
highest fluences, suggesting a dichotomy between the low 
and high fluence regimes. 

These arguments provide a framework for understand- 
ing the rest of the available phase diagram. At lower hole 
doping levels (i.e., x = 0.2, x = 0.3) and higher excitation 
fluences, the data represent the normal-state behavior 
seen at higher dopings. In the low fluence regime, how- 
ever, these samples do not exhibit the density-dependent 
recombination seen for the optimally and overdoped sam- 
ples, precluding an analysis to determine Fermi surface 
topology. Rather, these data are qualitatively explained 
by a simple model that accounts for the spatially inte- 
grated nature of the pump-probe measurement. 

In the Co-doped samples, low-fluence measurements 
below T c in an underdoped (y = 0.12) sample indicate 
quasiparticle recombination that is weakly dependent on 
incident fluence. This intensity dependent recombina- 
tion disappears as further electron doping submerges the 
innermost hole pocket below the Fermi level; for large 
values of Co-doping, the only observed relaxation is the 
slow decay we previously assigned to the outer hole band. 



II. EXPERIMENTAL SETUP 

We employed a Ti:sapphire oscillator system which 
produces 60 fs pulses at central wavelength 795 nm (pho- 
ton energy 1.5 eV) at an 80 MHz repetition rate. After 
passing through a prism compressor designed to account 
for the dispersion of the setup, the repetition rate of the 
laser was reduced to 1.6 MHz with an electro-optic pulse- 
picker in order to eliminate the effects of cumulative heat- 
ing of the sample. The beam was then split into two paths 
with a 9:1 ratio between the pump and probe. The hor- 
izontally polarized pump beam was chopped at 102 kHz 
by a photoelastic modulator/polarizer pair for lock- in de- 
tection and then passed through neutral density filters in 
order to tune the incident pump fluence from 44 /zJ/cm 2 
to 44 nJ/cm 2 . After attenuation, the pump beam was 
passed through a fast scanning delay stage which swept 
the pump delay by ^40 ps at a rate of 30 Hz before being 
focused to a 60 /im FWHM spot-size on the sample. 

After being split from the pump, the probe beam 
passed through a computer-controlled motorized delay 
stage and was rotated to vertical polarization before be- 
ing focused to a 60 ^m FWHM spot-size on the sample 
in the same location as the pump. The reflected probe 
beam was separated from the pump by spatial and po- 
larization filtering, then focused onto a photodiode. The 
photodiode output was low-pass filtered and measured 
by a lockin amplifier synched with the photoelastic mod- 
ulator. We used a fast time constant (30 /its) so that the 
lockin output could be input to an oscilloscope triggered 
to the fast scanning stage. We averaged several fast-scan 
stage sweeps to record each data trace. 

All samples were grown using the self-flux method in 
proportions given by their chemical formula. A variety 
of sample characterization methods were employed with 
the most rigorous testing reserved for the optimally hole- 
doped (x — 0.4) sample as follows: the induction-coupled 
plasma (ICP) technique was used to determine the Ba/K 
ratio, the uncertainty of which was determined to be less 
than 29p^. Energy-dispersive X-ray (EDX) spectroscopy 
was used to determine spatial variation in sample dop- 
ing homogeneity. A plot of the resulting EDX analy- 
sis pattern may be found in Figure lb of Ref. [M] In 
a recent EDX analysis on a different batch of samples, 
the maximum spatial variation was found to range from 
x = 0.39 to x = 0.44. Resistivity and magnetic suscepti- 
bility measurements have been independently performed 
and presented^; the sharpness of the transition was used 
as a gauge of sample purity. These results are shown in 
Figs, la and lb of Ref. US 

A number of other studies on samples from the 
same batch of optimally-hole doped compo unds have 
been performed, including ARPES 12 14 28 , Infrared 
spectroscopy^, and thermal conductivity^ further in- 
formation regarding sample purity and homogeneity may 
be found therein. We performed SQUID magnetome- 
try on all the samples used in this study and measured 
AT C 1 K at each sample's transition temperature. 
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FIG. 1. Schematic phase diagram of the Ba-122 pnic- 
tides for the hole-doped Bai_ a! K :E Fe2As2 and electron-doped 
BaFe2- M Co; ; As2 samples used in this study. The T c values 
are denoted by black squares and the Tsdw values are given 
by red circles. All T c values were measured by SQUID mag- 
netometry. The values for Tsdw on the electron-doped side 
of the phase diagram were determined from the same batch of 
samples 37 , while the values for Tsdw on the hole-doped side 
of the phase diagram are from the literature^. 



III. RESULTS AND DISCUSSION 
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FIG. 2. Fast, intensity-dependent relaxation of the normal- 
ized reflectivity transients (AR/ R) norrn in the superconduct- 
ing state, (a) The decay rate of (AR(t) / R) norm in the op- 
timally doped compound [x = 0.4) systematically decreases 
with decreasing intensity at 7 K (from bottom to top: $ = 
0.7 ^J/cm 2 (light green), 2.2, 4.4, 7.0, 8.9, 12.8 to 16.1 ^J/cm 2 
(black)). Raw traces from the overdoped hole-doped samples 
Bai_ I K a; Fe2 As2 for (b) x = 0.5 and (c) x = 0.6 at T = 5 K are 
qualitatively identical to those for the optimally doped sample 
in panel (a). In both of these plots, the absorbed fluences are, 
from bottom to top, $ = 0.7 ^J/cm 2 (blue), 1.4, 2.2, 4.4, 7.0, 
12.7, 20.1, and 36.6 ^J/cm 2 (black), (d) (AR/R) norm mea- 
sured at four different repetition rates are identical, verifying 
the absence of cumulative heating ($ = 37 /iJ/cm 2 , T=7 K). 



A schematic representation of the Ba-122 phase di- 
agram is presented in Fig. [I] The data points in- 
dicate the T c and Tsdw values for the samples used 
in this study. Specifically, we examined hole-doped 
Bai_a;K x Fe2As2 with x = 0.2 — 0.6 in increments of 
0.1, spanning the regime from underdoping to overdop- 
ing. On the electron-doped side of the phase diagram, 
we studied BaFe2_j / Co, / As2 samples at doping levels of 
y = 0.08 — 0.20 in increments of 0.04, which spans un- 
derdoped to optimally doped samples. The parent com- 
pound was also measured. 

We start by describing our measurements in 
Bai_a:K x Fe2As2 in the optimally (a: = 0.4) to overdoped 
(x = 0.6) regime at the low excitation densities. Shifting 
our attention to the high fluence regime provides a frame- 
work for measurements in the underdoped x = 0.2, 0.3 
samples, which we follow with a discussion of measure- 
ments in the parent compound and BaFe2_ K Co a; As2 at 
all dopings. 



A. Bai_ a; K :E Fe2 AS2: The Low Fluence Regime 

Characteristic short-time traces of the normalized 
change in reflectivity AR/R as a function of time at vari- 
ous fluences $ are shown in Fig. [2^, for the x = 0.4 sample 



at T=7 K. There is an initial decrease of the reflectivity 
at the arrival of the pump beam at time t = 0. Recovery 
to equilibrium depends strongly on the incident fluence 
with higher fluences exhibiting a faster initial relaxation 
rate than lower ones. Qualitatively identical intensity 
dependence is observed in the x = 0.5, 0.6 samples, as 
shown in Figs. [2]d and [2J;. In all cases, higher fluences 
are marked by faster initial relaxation of the signal than 
for lower fluences. We note the presence of the beginning 
of a heavily damped oscillation in the highest fluences of 
this Figure (not shown for the x = 0.4 sample, although 
it is was observed) due to Stimulated Brillouin Scatter- 
ing (SBS) which disappears with decreasing fluence. This 
characteristic of the data will be described in further de- 
tail below within the context of the high-fluence regime 
measurements. 

In order to rule out steady state heating as the source 
of this intensity dependence, we used a pulse picker to 
vary the repetition rate. Figure [2[l shows 7 K transients 
taken with an absorbed pump fluence $ = 37 /iJ/cm 2 at 
repetition rates ranging from 200 kHz to 1.6 MHz. We 
observe no discernible change in the recovery dynamics. 
Measurements performed above T c exhibited an absence 
of the intensity dependence of Figs. [2^,-c, verifying its 
origin in superconductivity^. 
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FIG. 3. Slow relaxation of (AR/R) norm in the supercon- 
ducting state, (a) (AR/ R) norm at 7 K and at room tem- 
perature. The slow, excitation density independent portion 
(black) switches off to become a nearly flat response in the 
normal state, which clearly shows the first oscillation due to 
stimulated Brillouin scattering, (b-d) A more complete view 
of the temperature dependence of (AR/ R) norm near T c ob- 
tained at $ = 4.4 J /cm 2 in the x = 0.4 — 0.6 samples. We 
note a sharp decrease in the offset across the transition (i.e. 
37 K to 42 K) due to switching off of the long-time component 
with loss of superconductivity. Above T c there is an upturn of 
the signal evident at short times due to stimulated Brillouin 
scattering^. 



The reflectivity transients of Figs.[2^,-c tend to an offset 
at the end of the 30 ps measurement time window. This 
offset is the beginning of a slow, intensity-independent 
decay with a characteristic decay time of ^500 pa^, rep- 
resented for the x = 0.4 sample at 7 K by the black 
trace in Fig. [3^,. As with the fast, intensity dependent 
component, the slowly decreasing offset is observed to 
abruptly shut off for all samples at their respective tran- 
sition temperatures, leading to a much longer decay of 
lower magnitude, as seen in the red trace of Fig. [3^. This 
switching off of the long-time component was also a uni- 
versal characteristic of the x = 0.4 — 0.6 samples, which 
is demonstrated by Figs.[3]3-d. 

Further insight into the data is obtained by consid- 
ering the initial amplitude of the reflectivity transients. 
Fig. ^presents (AR(0)/R) as a function of the absorbed 
pump fluence in the low-fiuence regime at T = 7 K in the 
x = 0.4—0.6 samples. Lines of slope 1 are shown with the 
data, indicating that at these fluences, AR(Q) /R oc uq is 
directly proportional to $. This linear proportionality 
signifies that these experiments represent the low-flucncc 
regime. Significantly, a linear dependence of uq on $ 
is an indication that the photoinduced excitations are 
fully gapped; the number of excitations is simply propor- 
tional to the laser energy absorbecP4 In the presence 
of a line node in the gap of the probed excitations, a 
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FIG. 4. (a) \(AR/R) max \ plotted as a function of <E> at 
T =7 K for the x — 0.4 — 0.6 samples, each offset for clarity. 
Dashed lines show a slope of one. Linear behavior in the low 
fluence regime is an indication that the optically accessed ex- 
citations are fully gapped, as described in the text, (b) The 
absolute value of the maximum change in reflectivity versus 
both $ and the average energy absorbed within a penetration 
depth for the x = 0.4 sample. The red curve is a fit to a sim- 
ple saturation model yielding $ sat = 25 /iJ/cm 2 at an initial 
temperature of 7 K. 



linear dependence of the density of states g{E) on en- 
ergy E produces an estimate for the total energy ab- 
sorbed by the quasiparticles as Et ~ / g(E)EdE ~ E 3 . 
In such a case, the total number of excited quasipar- 
ticles riph is given as n p h ~ / g(E)dE ~ E 2 . Thus, 
the number of quasiparticles in terms of the energy de- 

2/3 

posited is given as n p h ~ E T . Thus, a line node in the 
vicinity of the photoexcited quasiparticles should yield a 
sublinear dependence of the initial reflectivity transients 
AR(0)/R oc $ 2 / 3 , which is not observed. 

With increasing <f>, AR(0)/R displays saturation be- 
havior at high excitation densities, as seen for the x = 0.4 
sample 7 K data in Fig. |4Jd. We fit these data to a sim- 
ple saturation model which accounts for the exponential 
penetration of the light into the sample, given by 

AR(0)/AR sat = A" 1 / e-*l\l - e~^^)dz (1) 



where 3>(z) = $(0)e~ z / A is the laser fluence at a depth 
of z beneath the sample surface, AR sat is the change in 
reflectivity at saturation and A is the optical penetration 
depth (26 nm in optimally doped Bai_ x K x Fe2AsJ^). 
The fit yields a saturation fluence of $ sa t = 25 /iJ/cm 2 
and AR sat / R = 1-34 x 10~ 3 . Below, we restrict our 
analysis to $ < 2.5 /iJ/cm 2 which is clearly in the linear 
range of Fig. |4Jd for all three dopings. 

In the low-excitation regime, the dynamics of the pho- 
toinduced quasiparticles may be understood within the 
framework of the Rothwarf- Taylor model. This phe- 
nomenological model was originally developed to describe 
nonequilibrium quasiparticle recombination in tunneling 
experiments^. Here, the coupled recovery of a nonequi- 
librium concentration of quasiparticles n and their bind- 
ing boson N may be described by the rate equations 
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FIG. 5. Decay-rate analysis of data with a view towards gap 
symmetry, (a) The initial decay rate {'Jro) as a function of 
pump fluence at 7, 8 K for x — 0.4 — 0.6. The linear depen- 
dence of recombination rate on initial population is a hallmark 
of bimolecular recombination. Dashed lines show linear fits 
to data. The thermal decay rate at each temperature was 
obtain by extrapolating the fits to zero fluence. There is little 
significant difference between the various dopings. (b) The 
thermal relaxation rate 7t/i plotted as a function of tempera- 
ture for x = 0.4—0.6. The dashed line shows ~ T 2 dependence 
as a guide along for all doping levels, suggesting a negligible 
change in Fermi surface topology between differently doped 
samples. 



N = I ph + Bn 2 /2 - lpc N -(N- N eq ) le 



(3) 



where, I qp and I p h are the external quasiparticle and bo- 
son generation rates, respectively, j pc is the pair creation 
rate via annihilation of gap energy bosons, N eq is the 
equilibrium boson number density, 7 esc is the boson es- 
cape rate, and B is the bi-molecular recombination con- 
stant. The rates Bn, 7 pc and 7 esc describe the three dif- 
ferent relaxation pathways for the energy deposited into 
the system by the ultrashort pulse. 

Depending on the relative magnitudes of these three 
rates, the solutions of these equations display different 
characteristics. In the simplest case, if 7 pc <C Bn or 
7esc 3> 7pc then the pair creation term in the first 
equation can be ignored and the two equations decou- 
ple. In this limit, quasiparticles display simple bimolecu- 
lar kinetics governed by the recombination coefficient B, 
as is observed, e.g., in Figs. [2^,-c. In the limit where 
lesc ^ "/pa Bn, quasiparticles and bosons come to a 
quasiequilibrium and the combined population decays 
with a slow rate proportional to 7 esc . This regime, known 
as the phonon bottleneck, is consistent with the slow 
component of the signal, shown in Fig. [3^. 

In the decoupled regime, the RT equations permit de- 
termination of the thermal decay rate of quasiparticles 
('Jth) an d the recombination coefficient IPh This may be 
seen by explicitly drawing a distinction between photoin- 
duced quasiparticles (n p h) and thermally present ones 
(nth) and rewriting n = n p h + nth- Eq. [5] then becomes 



dn/dt = -Bn 2 h - 2Bn ph n th , 



(4) 



where the second term Bn p h,nth represents the recom- 
bination of photoinduced quasiparticles with their ther- 
mally populated counterparts. We consider the initial 



recombination rate 7 r o, defined as 

Tro = -0-/nph)(dn p h/db)\t-x) (5) 
= Bn ph + 2Bn th = l P h + Zjth (6) 

where ~fph,th = Bn p h,th are the photoinduced and ther- 
mal decay rates. 

Experimentally, we deduce j r Q from the initial slope 
of the reflectivity transients following the peak. As sug- 
gested by Eq. [6j we obtain 7^ by extrapolating the linear 
fits to zero excitation density (i.e., n p h oc $ — > 0) for each 
temperature, as shown in Fig. [5^, at 7 and 8 K. The slope 
of the fits is proportional to B, which we observed did not 
strongly depend on the temperature for any of these sam- 
ples. The intercept jth wa s observed to increase with T 
due to the greater thermal quasiparticle population. We 
shall examine each quantity in turn. 

In order to determine the value of B, an estimate of n p h 
for a given $ is required, which we will conveniently pick 
as the saturation fluence. The average energy density de- 
posited by the pump beam within the penetration depth 
at saturation is $ sa t/A = 9.6 J/cm 3 , using the value for 
$ sat determined above. We compare this value with an 
estimate based on the assumption that saturation occurs 
when the quasiparticle system receives an energy compa- 
rable to the condensation energy, which we approximate 
from the BCS theory as l/2N(E f )A 2 . Here, N(E f ) is 
the density of states at the Fermi level and A is the su- 
perconducting gap energy. Using an average gap size of 
A=10 meV and N(E f ) = 7.3 eV _1 per unit celP 1 , we 
obtain a condensation energy of 0.3 J/cm 3 . This value 
is significantly smaller than the experimentally measured 
saturation energy of 9.6 J/cm 3 . Hence, only 3% of the 
energy goes into creating quasiparticle pairs. Using this 
value to convert the laser fluence into excitation density, 
we obtain B = 2.48 x 10~ 9 cm 3 /s. For a BCS supercon- 
ductor in the dirty limit, B is proportional to the ratio of 
the electron phonon coupling constant to the density of 
states at the Fermi energy^. In the case of pnictide su- 
perconductors with multiple gaps, a detailed theoretical 
model for the quasiparticle recombination is necessary to 
further interpret the measured value of B. 

Figure [5]d presents jth as a function of temperature for 
x = 0.4 — 0.6. Significantly, the data presented in this 
Figure reveal that the thermal decay rate 7^ = Bnth 
is proportional to T 2 below T c for these doping levels. 
This suggests the presence of a node in the gap around 
which nth ~ T 2 , as opposed to the exponential depen- 
dence expected from an isotropic gap. However, the lin- 
ear dependence of the magnitude of AR/R on $ at weak 
excitation levels (Fig. [4^,) argues against the presence of 
nodal quasiparticles, which would lead to AR/R oc $ 2 / 3 , 
as discussed above. 

The phenomenological consistency of the pump-probe 
measurements in the doping range x — 0.4 — 0.6 indi- 
cate that the physics of these materials is consistent with 
the band-dependent picture we have developed for the 
x = 0.4 material. As this picture is central to our un- 
derstanding of nonequilibrium quasiparticle dynamics in 
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FIG. 6. Schematic representation of the pnictide band struc- 
ture. While five bands have been observed by ARPES, the 
inner hole band is nearly doubly degenerate and so is repre- 
sented by the single band a. 



the pnictides, we briefly summarize it using the simplified 
representation of a pnictide band structure in Fig. [6j 

LDA calculations^ have shown a lack of states 1.5 eV 
above and (accounting for renormalizatiorfSJ) below the 
Fermi level at the M point of the Brillouin zone, implying 
that the probe wavelength used in these experiments is 
not able to couple to quasiparticle dynamics at this loca- 
tion in the zone. This point is bolstered by an estimate 
of the Sommcrfcld parameter, which shows the density of 
states of the hole bands to be several times larger than at 
the electron pockets^, resulting in a much greater num- 
ber of photoinduced quasiparticles at the hole pockets as 
compared with the electron pockets. 

Thus, the electron bands 7, 5 do not contribute signifi- 
cantly to signal, i.e., AR(t)/R derives almost exclusively 
from the hole bands. In particular, strong couplin g be- 
tween the carriers in the a bands with those in 7, JHUll 
leads to the observed excitation density dependent dy- 
namics, while the isolation of carriers in (3 from the other 
bands is represented by their slow, bottlenecked decay. 
Furthermore, since the entirety of the signal is due to 
the hole bands, we may thus interpret the observation 
AR(Q)/R (x $ to indicate that it is these bands which 
possess fully gapped excitations. 

As the excitation density is lowered, the thermal quasi- 
particle population becomes comparable to, and then ex- 
ceeds, the photoinduced one. Nonequilibrium quasiparti- 
cle recombination is then dominated by the second term 
in Eq. [4j which represents a photoinduced quasiparticle 
recombining with a thermal one. This measurement thus 
becomes particularly sensitive to the presence of either a 
highly anisotropic gap or a node on the Fermi surface, as 
the locally small gap energy leads to a relatively higher 
concentration of thermal quasiparticles in its vicinity. 

Here, the observation that jth ^ T 2 may be indicative 
of either a highly anisotropic gap or node on the elec- 
tron pockets, as interband scattering may be dominated 
by photoexcited carriers in a recombining with a rela- 
tively large population of thermally present carriers in 
7, 5. This scenario is compatible with an s± order pa- 
rameter with impurity scattering to account for the T 2 
behavior, although it does not require it. The fact that 
jth T 2 does not change up to doping levels of x = 0.6 



-1.0 



8 K, x=0.4 




50 100 
Time (ps) 



200 400 

Time (ps) 



600 



FIG. 7. Dynamics of the decay measured for a longer time 
window at T = 8 K and $ = 37 fj, J /cm 2 . Highly damped 
oscillations are observed with a period of ~ 40 ps (inset). 



suggests that the Fermi surface topology is consistent up 
to this doping. 



B. Bai„ I K I Fe2As2: The High Fluence Regime 

We now shift our attention to phenomena that require 
a different framework than that developed above, i.e., the 
observation of coherent oscillations in all samples at the 
highest fluences ($ > 10/iJ/cm 2 ) and the dynamics of 
the underdoped hole-doped samples where x = 0.2, 0.3. 

At fluences greater than 10/^J/cm 2 , we observe the on- 
set of highly damped oscillations in the short time tran- 
sients shown in Figs[2]3 and [2]: that become more appar- 
ent with increasing fluence. A clearer depiction of this 
effect is provided by Fig. [7] by extending the measure- 
ment window for the x — 0.4 sample. These oscillations 
are due to stimulated Brillouin Scattering^, where inter- 
ference between the portion of the probe beam reflected 
from the sample surface and from the propagating strain 
pulse launched by the pump modulates the signal in the 
time domain. The observed period of ~ 40 ps is consis- 
tent with the expected value of X pr / (2nv s cos 0}pD where 
A pr is the wavelength of the probe, 9 the angle of inci- 
dence, v s the speed of souncP^and n the refractive index. 
The fast damping rate T is set by the small optical pen- 
etration depth of the light (r = v s /X). 

The strength of the SBS signal is linearly proportional 
to the thermal expansion coefficient. Thus, the disap- 
pearance of the oscillations in the low fluence data below 
T c is a manifestation of its large drop at the transition 
temperature T c for all Ba-122 pnictides. This drop is due 
to the thermodynamic relationship between c p (T) and 
the expansivit y 3 1 " 33 l and is more pronounced for expan- 
sion along the c-axis as a result of the layered nature of 
the pnictides^. The presence of an acoustic wave may be 
an indication that the sample has been partially driven 
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FIG. 8. Raw traces in the underdoped hole-doped samples 
Ba i _ :r K ;c Fe2As2 for (a) x = 0.2 (T c = 22 K, T S dw ~ 100 K) 
and (b) x = 0.3 (T c = 22 K) at T = 5 K. The absorbed 
fluences are $ = 36.6 (black - (a) only), 20.1, 12.7, 7.0, 4.4, 
2.2, 1.4, and 0.7 (blue) /uJ/cm 2 . Also shown are simulated 
results for AR(t) /R at various values of $ using a simple two- 
component model in the (c) x — 0.2 sample and the (d) x = 
0.3 sample at the same fluences. The simulations, which are 
performed without free parameters, show excellent qualitative 
agreement with the data. 



into the normal state, even at the relatively modest flu- 
ences used here (< 37 /iJ/cm 2 ). 

We examine the plausibility of this scenario by simple 
energy conservation arguments. We recall that the exper- 
imentally measured saturation energy, estimated above 
as 9.6 J/cm 3 , supports the hypothesis that high fluences 
deplete the superconducting condensate. However, this 
also indicates that a large fraction of the energy deposited 
by the laser is transferred into other internal degrees of 
freedom or is carried out of the probed region by through- 
plane transport of the excited electrons. The large exci- 
tation gradient that results from the short penetration 
depth of the pump may result in fast diffusion of the ex- 
cited electrons, quickly transporting a significant fraction 
of the absorbed energy out of the experimentally probed 
region. 

Shifting focus to the underdoped samples, Figs. [HJi 
and [8]d show the T = 5 K reflectivity transients in the 
x = 0.2 and x — 0.3 samples. In the x = 0.2 material, 
T c = 22 K and a spin-density wave transition occurs at 
Tsdw = 120 K. In the x = 0.3 sample, Tsdw is sup- 
pressed to 100 K and T c is elevated to 32 KP. The dy- 
namics observed here are markedly different than those 
of the optimally and overdoped samples. In the x = 0.2 
compound, the high fluence traces resemble the normal 
state of the x = 0.4 — 0.6 samples. With reduced flu- 
ence, we see the emergence of an intensity-independent, 
slowly-decaying positive component reminiscent of the 
slow dynamics of Fig. [3^,. 



The dynamics of the x = 0.3 sample shown in Fig. [8] 
reveal behavior similar to those of the x = 0.2 sample 
at the highest fluences, again in contrast with the mea- 
surements of Figs. [2£i-c. Here, we observe that with a 
reduction of the incident laser fluence, a positive compo- 
nent of the signal emerges whose recovery dynamics are 
intensity independent and decay on the relatively quick 
~ 10 ps timescale. 

The absence of a bimolecular intensity dependence may 
be interpreted within the framework of the Rothwarf- 
Taylor equations as bottlenecked recombination. Unfor- 
tunately, this precludes an analysis to determine gap 
symmetry. However, we may understand these data 
by applying a simple model which considers the quasi- 
equilibrium temperature of the sample T(z) as a func- 
tion of depth in the sample z after excitation, the results 
of which are shown in Figs. [8b and [§ji. Conservation of 
energy provides an expression for T(z): 



— exp (—z/d) 
a 



T{z) 



c p (T)dT. 



(7) 



To 



Here, the specific heat c p (T) was modelled as a fifth- 
order polynomial using published data on the x = 0.4 
compound^. Equation [7] was then integrated and in- 
verted numerically to reproduce the initial temperature 
profile, depicted in Fig. [9^,. We note that high fluences 
produce a significant increase in temperature. In con- 
trast, this increase is very modest at the low fluences 
examined in Sec. MI Al 

The pump-probe measurement is an integrated mea- 
surement in depth. We therefore calculate the fraction A 
of the sample which has been brought above the transi- 
tion temperature from 



A 



f™Q(z')exp(-2z'/d)dz' 
/ oo exp(-2z'/d)dz' : 



(8) 



where Q(z') = 1 if T(z') > T c and = otherwise. The 
factor of two in the above equation accounts for the reflec- 
tion of the beam back through the sample. In practice, 
this integral was numerically evaluated to 3 penetration 
depths. The result is shown in Fig. [9]b. For an initial tem- 
perature of 7 K, the entire sample remains below T c for 
$ < 6 /iJ/cm 2 and achieves a 90:10 distribution between 
normal and superconducting state at $ < 19.8 /iJ/cm 2 . 
The small discrepancy between the value obtained from 
this simple model and the experimentally measured sat- 
uration fluence may be due to through-plane electronic 
diffusion mentioned above. 

We consider the signal to comprise two parts. The 
first is a fast decaying negative component of timescale 
Tf with an offset c as f(t) = —(exp (— t/rj) +c) to repre- 
sent the normal state contribution. The second is a slow, 
positive component for the superconducting contribution 
of timescale t s as s(t) — exp(— t/r s ). The two compo- 
nents are then added in the proportion Af(t) + (l— A)s(t), 
where the coefficient A is defined as the fraction of the 
sample seen by the probe beam with T(z) > T c . Here, 
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FIG. 9. (a) Simulated temperature profile as a function of 
depth for the fluences used in this study. The dotted lines 
indicate the values of transition temperature for the x = 0.2 
(T c = 22 K) and the x = 0.3 (T c = 32 K) samples, (b) 
Fraction of the sample heated into the normal state as inter- 
rogated by the penetration depth of the probe for the x — 0.3 
sample. 



we compute A as a weighted average of the temperature, 
as in Eq. [8] 

The results of this calculation using c = 0.2, 77 = 1 ps 
and t s = 400 ps are shown in Fig. [HJ; for the x = 0.2 sam- 
ple. Fig. [HJl uses the same vales for c and 17 but a value 
of t s = 10 ps for the x = 0.3. In both cases, we observe 
excellent qualitative agreement between the experimen- 
tal data and the model using no fitted parameters. Data 
acquired at higher initial temperatures showed similar 
agreement with this model for both x = 0.2 and x = 0.3. 

We note that this model is not able to replicate the 
dynamics of the optimally doped and overdoped samples 
discussed in the previous section. Even when account- 
ing for the change in sign of the phenomenological su- 
perconducting component s(t) from positive to negative, 
no combination of simulation parameters can produce an 
initial decay rate that depends linearly on excitation flu- 
ence. This provides further evidence that the intensity 
dependence observed in the x = 0.4 — 0.6 samples arises 
from bimolecular recombination kinetics. We thus em- 
phasize that this analysis has no impact on the analy- 
sis performed in the low-fluence regime. First, in those 
traces, we deduced the gap symmetry from the initial 
quasiparticle decay rate 70, which was obtained from the 
slope of the reflectivity transients immediately after the 
peak, i.e., t < 1 ps. This is much shorter than the time 
required for an equilibrium temperature to be established 
between the electronic system and the phonons. Second, 
all the data used in the analysis of Sec. |III A| were at 
such low fluences that the simulated equilibrium sample 
temperature did not depart significantly from the initial 
value. 

The physical processes at play in the x — 0.2 and 
x = 0.3 samples may also be understood within the 
framework of the band-dependent dynamics discussed 
in the previous section. Ultrafast measurements on the 
Bai_ 2; K a ;Fe2As2 samples for values of x < 0.3 indi- 
cate a competition for itinerant electrons between the 
SDW state and the superconducting state^ 5 . This ob- 
servation is consistent with ARPES measurements on 
Bao.75Ko.25Fe2As2 which have shown that the supercon- 
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FIG. 10. Raw traces in the underdoped hole-doped pnic- 
tide Bao.7Ko.3Fe2As2 (T c = 32 K) above T c . In (a), we note 
that there is some residual intensity dependence of the signal, 
where the initial recovery of the fast component is quicker at 
the highest fluences than at lower ones. This is identical to 
the superconducting state, despite the fact that at T = 34 K, 
the sample is in the normal state and the SDW state does not 
exhibit intensity dependent recombination (see below). After 
the initial electronic recovery, the amplitude of the acoustic 
wave is diminished at lower fluences. (b) In contrast, the 
fluence dependence of the electronic recovery has all but dis- 
appeared at 60 K. In both plots, the absorbed fluences are <3? 
= 36.6, 20.1, 12.7, 7.0, 4.4, 2.2, 1.4, 0.7, and 0.44 /ij/cm 2 . 



ducting fraction in a and 7, S are suppressed due to their 
contribution to SDW order via Fermi surface nesting^. 
At x = 0.2, this nesting is even slightly more favored and 
hence the bulk of the optical pump-probe response in the 
superconducting state arises from the slow, intensity in- 
dependent decay in the outer hole band /3 with little to 
no contribution to the signal from the other bands. 

In the x = 0.3 sample, the absence of intensity de- 
pendence is again due to the presence of a SDW state 
between a and 7, 5 which competes with superconduc- 
tivity for carriers in these bands. It is, however, unclear 
at present why the dynamics change from slow to fast at 
this particular doping. 

The presence of a pseudogap state in the 1111 
system has already been argued for based on ultra- 
fast pump-probe measurements^ ARPEyHI and optical 
conductivity^ have indicated the same for underdoped 
Bai_ 2; K a .Fe2As2 samples. In the context of the intensity 
dependent dynamics studied here, pseudogap behavior is 
suggested by the intensity dependence present above T c 
in Fig. 3c of Ref. [27] for x = 0.4. A clearer demonstra- 
tion of such fluence dependence is depicted in the time- 
resolved reflectivity transients of Fig. [TOfc ,, which were 
measured in the x — 0.3 sample 2 K above the transition 
temperature. There is a discernable tendency towards 
faster initial relaxation of the reflectivity transients with 
an increase in temperature. This disappears by 60 K, as 



may be seen in Fig. 10 :>. 

In light of the band-dependent relaxation dynamics 
discussed above, we posit that those carriers participat- 
ing in the pseudogap state must derive from the a band 
because these carriers give rise to intensity dependent re- 
combination in the superconducting state. This conclu- 
sion is consistent with prior ARPES measurements on 
Bao.75Ko.25Fe2As2 which show the presence of a pseudo- 
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FIG. 11. (a) Schematic representation of the band struc- 
ture of the Ba-122 pnictides. Although this is somewhat of a 
simplification that ignores the change in k z dependence with 
doping^, the predominant effect is a shift of the Fermi level 
up with the addition of electronp^ in the parent compound, 
a and 7, S nest, leading to SDW behavior. In the underdoped 
regime (y = 0.12), a does not nest appreciably with 7, S. With 
the addition of more holes, the a band is submerged below 
the Fermi level when the doping reaches y — 0.15, as mea- 
sured by Sekiba et al.™ (b) Raw (AR/ R) norm traces for the 
parent compound at 20 K. The moderate amount of intensity 
dependence in these traces at short times is due to acoustic 
suppression from the diminished thermal expansion parame- 
ter at lower temperatures/fluences. (c) Reflectivity transients 
in the underdoped sample (y = 0.12, T c = 12 K) show a mild 
amount of intensity dependence of the faster relaxation with 
increasing fluence, similar to that seen in the hole-doped sam- 
ples, (d) Fluence dependence of the optimally doped sample 
(y — 0.2, T c — 22.5 K) indicates entirely different relaxation 
dynamics with fluence. While there appears to be intensity 
dependence present here, the slope of the initial transient, 
70 is identical for all fluences. The differences between the 
various fluences and the origin of this non-uniformity of the 
relaxation dynamics are due to a temperature effect, as de- 
scribed in the text. 



gap in a but not in the other bands^. 



BaFe 



2-1/ 



Co H As2 



We now shift our focus to the electron doped side of 
the phase diagram. Fig.fTTfc, depicts the qualitative effect 
of doping on the Fermi level. As electrons are added to 
the system, the Fermi level rises in a rigid band manner, 
as has been observed by Brouet et alPS At a Co-doping 
of y = 0.15, the a band begins to submerge below the 
Fermi level at various points along fcj^. By the point of 
optimal Co-doping at y = 0.20, a has been submerged 
well below the Fermi level. 

The reflectivity transients of Fig. [TT] exhibit the con- 
sequences of the disappearance of a. As a point of com- 



parison, Fig. |llp shows data from the parent compound 
taken at 20 K, well below the SDW transition temper- 
ature Tsdw = 143 K. None of the intensity dependent 
relaxation dynamics associated with the x > 0.4 hole- 
doped samples were observed in the undoped sample. We 
note the small effect of the suppression of the thermal 
expansion predominantly along the c-axis^^U which in- 
hibits the generation of an acoustic wave. Qualitatively 
identical behavior was observed in an underdoped sam- 
ple (y = 0.08, not shown) for which the SDW transition 
is suppressed to 95 K but which does not superconduct 
at low temperatures. 

Fig. [TT] ; shows raw data for BaFei.8sCoo.i2As2 
(T c = 12 K) taken at different fluences ranging from 
0.14 /iJ/cm 2 to 13 ^J/cm 2 . As with the x = 0.4 - 0.6 
samples, we observe a fast intensity dependence of the 
signal at short times with higher laser intensities corre- 
sponding to faster initial recombination rates. We note 
that the relatively low T c — 12 K of the y = 0.08 sample 
does not allow us to use the initial decay-rate analysis 
employed above to determine Fermi surface topology. 

With the submergence of the inner hole band a at 
y = O.lfP^, the characteristics of (AR/R) norm change 
markedly. Although qualitatively similar behavior was 
observed in the y = 0.16 (T c — 16 K) sample, these 
dynamics are most clearly represented by the y — 0.20 
sample due to its slightly higher T c — 22 K. As may be 
seen in Fig. [Tip ., AR/ R traces scaled to the initial elec- 
tronic response at t = are qualitatively identical (within 
a sign) to those of the x = 0.2 underdoped hole-doped 
sample in Fig[8ju Indeed, we observe that the polarity of 
the signal has changed from the underdoped to the opti- 
mally doped side of the phase diagram. A similar effect 
has also been observed in the Bi 2 Sr2Cai-vDyyCu2 08+d 
(BSCCO) system of cuprate superconductors^. 

In order to better interpret the dynamics of Fig. [Tip 1 , 
we modelled the heat flow out of the excited region after 
the establishment of an initial equilibrium temperature. 
This was done by numerically solving the heat diffusion 
equation using literature temperature-dependent thermal 
conductivity^ and heat capacitjJ^data while assuming a 
temperature-independent density. The initial condition 
was specified using Eq. [7] and the incident fluence. In 
the limit of bottlenecked recombination, AR(t)/R should 
be independent of initial quasiparticle density once the 
entire probed region is below T c . This indicates that 
the data should be scaled at the time t when the entire 
sample becomes superconducting. 

Reflectivity transients were acquired out to 600 ps and 
then rescaled at the point in time at which the entire 
probed depth is below T c . This rescaling is shown in 
Fig. |12^ ,, where we observe a remarkable convergence 
of the data at long times. The highest fluences mea- 
sured produce the normal state response, i.e., an initial 
"spike" followed by a quasi-steady state that does not 
decay within the time window. As the excitation fluence 
is reduced, the initial recovery gradually transitions to 
bottlenecked quasiparticle recombination at earlier and 
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FIG. 12. Evolution of (AR/R) norm in BaFe1.sCoo.2As2 for 
absorbed fluences of $ = 36.6 (black), 20.1, 12.7, 7.0, 4.4, 2.2, 
1.4, 0.7, and 0.44 (blue) /ij/cm 2 . As the fluence is increased, 
the recovery to equilibrium seems to slow due to a greater 
fraction of the probed sample being heated above T c , espe- 
cially at these long time scales. Proper scaling of the various 
traces yields a common, bottlenecked decay, reached by the 
sample when the entire probed region has cooled below T c 
due to thermal diffusion, (b) Comparison of long-time traces 
of the optimally-doped hole-doped sample (Bao.6Ko.4Fe2As2) 
and the optimally electron-doped (BaFe1.gCoo.2As2) material 
at fluences of 0.7 and 0.4 /nJ/cm 2 . Even though the Fermi 
surface topologies and gap energies of the two samples differ 
markedly, we note that the relaxation dynamics are nearly 
identical, indicating a common underlying relaxation mech- 
anism, i.e., bottlenecked quasiparticle recombination in the 
outer hole band (/?). This slow relaxation switches off at the 
respective values of T c of both samples. 



earlier times. At the lowest fluences (<& = 0.7 /iJ/cm 2 
and <f> = 0.44 /zJ/cm 2 ), the entire probed region is below 
T c and there is an intensity independent, bottlenecked 
relaxation of the signal. 

This behavior recalls the slow dynamics of the opti- 
mally hole-doped sample shown in Fig. [3^. A comparison 
of these two data sets is provided in Fig. [T2p which plots 
data from the optimally hole-doped sample (x — 0.4) at 
<f> = 2.2 /iJ/cm 2 alongside data recorded at the lowest flu- 
ences in the optimally electron-doped sample (y = 0.2). 
We observe nearly identical behavior, suggesting that 
these decays arise from the same set of dynamics. As 
mentioned above, this slow decay was observed to switch 
off in both samples abruptly at their respective supercon- 
ducting transition temperatures, indicating that it origi- 
nates from superconductivity. 

We posit that the similarities between the long-time 
reflectivity transients produced by Bao.6Ko.4Fe2As2 and 
BaFe1.gCoo.2As2 may also account for the differences be- 
tween the data in Figs. [TT| ; andpTp!. That is, the disap- 
pearance of bimolecular recombination with the doping 
of electrons derives from the band-dependent recombina- 



tion dynamics presented above in Sec. Ill A The addition 
of electrons disrupts the near-perfect nesting of a with 
7, 5, weakening the effect of the magnetic resonance on 
Cooper pairing and resulting in a diminished intensity 
dependence of the y — 0.12 sample of Fig. [Tip . Further 
doping of electrons leads to the submergence of a and 
thus the remaining relaxation in the y = 0.20 sample is 
the slow decay of photoinduced carriers in /3. 



IV. SUMMARY AND CONCLUSION 

We have presented an ultrafast pump-probe study of 
the Ba-122 pnictides as a function of fluence and dop- 
ing on both sides of the phase diagram that provides 
evidence for band-dependent dynamics. For optimally 
doped Bai_ K K a ;Fe2As2, the initial recovery rate of the 
reflectivity transients depends linearly on the incident 
laser fluence $ in the low excitation limit. This com- 
ponent of the signal was argued to derive from the inner 
hole band a. A Rothwarf- Taylor analysis of this intensity 
dependence in the low-fluence limit yields a T 2 behav- 
ior for the thermal population of quasiparticles. At long 
times, AR/R tended towards a slow, intensity indepen- 
dent decay which originates from the outer hole band, 
f3. These observations were unchanged with an increase 
of K-doping up to levels of x = 0.6. Our results thus 
indicate fully gapped hole bands with nodal or strongly 
anisotropic electron bands in these samples. 

In underdoped samples (x — 0.2, 0.3), the intensity de- 
pendence disappears as electrons in a are taken up by 
SDW ordering and hence unable to participate in super- 
conductivity. The fluence dependence of AR(t) /Ris then 
governed by the decay in /3 in a proportion which may 
be estimated by considering the fraction of the probed 
depth which is driven above T c by the pump. 

In the BaFe2-yCoyAs2 samples, we observed intensity 
dependent recombination at low doping levels (y = 0.12) 
which completely disappeared with the submergence of 
the inner hole band a by the doping of electrons. The 
long-time dynamics of these samples were consistent with 
bottlenecked behavior due to quasiparticle recombination 
in the outer hole band (3. 

Intensity dependence was observed to persist above the 
transition temperature T c for samples with K-doping at 
levels of x < 0.4 and was taken as evidence for a pseudo- 
gap state. This pseudogap behavior was posited to arise 
from a, consistent with prior ARPES measurements. No 
such behavior is observed on the electron doped side of 
the phase diagram. 

Stimulated Brillouin scattering was observed in all 
samples above the superconducting transition tempera- 
ture but disappeared immediately below the transition 
temperature due to a large decrease of the thermal ex- 
pansion coefficient at T c . The presence of Brillouin Scat- 
tering below the transition temperature for the highest 
incident fluences used in this study (e.g., ~ 37 /iJ/cm 2 ) 
is an indication that, even at these relatively modest 
fluences, the optical pulse was able to bring the sam- 
ple into the normal state within the penetration depth 
of the excitation. At the highest fluences, this picture 
correlated well with basic modelling of reflectivity tran- 
sients in terms of a linear combination of a normal-state 
response and a superconducting, bottlenecked quasipar- 
ticle recombination, the relative proportions of which are 
dictated by a simple equilibrium temperature model. In 
the low fluence regime, however, where no part of the 
probed region could reach a temperature greater than 
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T c , either only bottlenecked or bare bimolecular recom- 
bination were observed. 
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